Cosmic Problems for Condensed Matter 
Experiment 



Tanmay Vachaspati 

CERCA, Physics Department, 
Case Western Reserve University, 
Cleveland, OH 44120, USA. 



Condensed matter analogs of the cosmological environment have raised the 
hope that laboratory experiments can be done to test theoretical ideas in cos- 
mology. I will describe Unruh's sonic analog of a black hole ("dumbhole") 
that can be used to test Hawking radiation, and some recent proposals on 
how one might be able to create a dumbhole in the lab. In this context, I 
also discuss an experiment already done on the Helium-3 AB system by the 
Lancaster group. 

Cosmology is a unique science where observations of the current or re- 
cent universe are used to infer about the very early universe. Furthermore, 
the physical ideas used to extrapolate back to the early universe are based 
on ideas that were developed in a non-cosmological setting. Surprisingly 
most of these ideas work extremely well. Yet there are also ideas that have 
no experimental verification but are used on the basis of theoretical extrap- 
olation. "Cosmology in the lab" is an emerging research area to search for 
analogs of the cosmological environment to enable experimental tests of ideas 
in cosmology. 

Cosmological problems are of two types. The first type contains prob- 
lems that are based on physics in which the cosmological environment is not 
crucial. It is just that the particular problem is more relevant to cosmology 
than to some terrestrial system. The outcomes of phase transitions fall into 
this class of problems. In cosmology, remnants of a phase transition, such as 
topological defects, are crucial since they could be observed today and lead 
to important clues about an earlier epoch. 

Problems of the second type are strongly based on the gravitational 
environment. For example, these could involve quantum effects in the pres- 
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ence of a "horizon" . In models of cosmic inflation, quantum field theoretic 
effects on scales larger than the cosmic horizon are supposed to have led 
to density fluctuations that then grew into galaxies. However, these quan- 
tum field theoretic effects are built upon the success quantum field theory 
(QFT) has had on very small (atomic) scales. Systems with horizons have 
not been constructed in the lab so far, and there are no experimental tests 
of QFT on superhorizon scales. Another system with a horizon is a black 
hole. Particles within the black hole cannot classically escape to the region 
outside the horizon. However, quantum effects are supposed to modify this 
picture, leading to Hawking radiation from the black hole 1 . It would be 
very desirable to find a way to experimentally test these ideas. 

Research on "cosmology in the lab" should be understood strictly as 
an attempt to simulate the cosmological environment in laboratory systems. 
It should not become an attempt to equate cosmological (or gravitational) 
phenomenon to some feature of a particular laboratory system. For example, 
cosmology in the lab cannot say if there is an underlying atomistic struc- 
ture to spacetime. However, if there was an atomistic theory of spacetime, 
condensed matter experiments might prove to be useful analogs. From my 
perspective, cosmology in the lab is similar to numerical studies of systems 
carried out on a computer. If the simulation is accurate, experiments can tell 
us things about cosmology that we would otherwise have no hope of finding 
out. 

In 1981 Unruh^ proposed an experimental analog of a black hole. He 
considered a fluid that is flowing at subsonic speed in the upstream region 
and at supersonic speed in the downstream region (see Fig. [J). Any fish 
in the downstream region cannot send sound signals to the fish in the up- 
stream region. Therefore there is a sonic event horizon. Classically, the fish 
upstream will receive no sound from the sonic hole i.e. the hole is "dumb". 
With quantum effects taken into account, however, there will be Hawking 
radiation of sound from the horizon of the dumbhole just as there is Hawking 
radiation of light from a black hole. 



To put Unruh's dumbhole on a more quantitative footing, the fluid is 
assumed to obey the irrotational, inviscid, Navier-Stokes equations: 

V x t> = , ^ + V • (pv) = 
at 

dv 1 

h v ■ Vv = — Vp- V$ (1) 

at p 



where denotes an external potential (e.g. gravitational potential). The 
pressure (p) is assumed to be a function of the density (p) alone. Fluctuations 
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subsonic flow 




Fig. 1. Unruh's vision of a waterfall as a sonic black hole or "dumbhole". The fish 
in the subsonic flow region cannot hear the screams of the fish in the supersonic 
region because the emitted sound travels too slowly to propagate upstream. The 
fish upstream sees a sonic horizon at the location where the fluid velocity becomes 
supersonic. 



around some background flow will correspond to sound. If po an d vq denote 
the background flow, 

p = po + bp 



v = v + Vc, 



Then it can be shown that 



1 



-9 







where the metric, g^ u , experienced by the sound is given by: 
ds 2 = (c 2 - vl)dt 2 + 2v dtdr - dr 2 - r 2 dtt 2 



(2) 



(3) 



This is the Painleve-Gullstrand-Lemaitre form of a black hole metric if we 
assume a spherically symmetric, stationary, convergent, background flow. 
The horizon is at vq = c s , that is, at the location where the fluid velocity 
equals the sound speed. 

The next step is to include quantum field theoretic effects. For this 
we have to go back to Eq. (J2J), find the mode functions, and compare their 
behavior near the horizon to that in the asymptotic region. This is the usual 
Hawking calculation and we will not show it here. We simply state Unruh's 
result for the Hawking temperature^ 



1 dv t 
2-7T dr 
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(4) 



where R is the distance over which c s changes. 
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Unruh's estimate highlights the difficulty of the problem. We need to 
accelerate the fluid by 300 m/s over a distance of 1 mm to get a temper- 
ature of a mere ~ 10~ 7 K. In units of the acceleration due to gravity, the 
required fluid acceleration is ~ 10 7 g! Achieving such an enormous acceler- 
ation is impractical given that we also want the fluid to remain very cold 
(/iK), so that we are able to detect the Hawking radiation. There are all 
sorts of instabilities that set in under such extreme conditions for all known 
systems ™. 

What is the alternative? A ray of hope emerges when one considers 
Visser's generalization B3 of Unruh's result. Unruh in his pioneering contri- 
bution only considered a fluid in which the sound speed is fixed. Visser 
generalized the setting to include the possibility that the sound speed may 
vary within the fluid. Then he found: 

r «*-(sb)!<*-">L <5> 

where c s = c s (t, x) is the sound speed and the derivative is evaluated at the 
location of the horizon (v = c s ). 

The ray of hope^is that we need not have huge gradients in v. Instead 
we might be able to arrange c s to vary very rapidly. An extreme case would 
be where v ~ everywhere but the properties of the fluid change at a bound- 
ary that is moving. The fluid is stationary with respect to the container but 
is moving with respect to the boundary. 

This was also the underlying idea in a proposal made by Jacobson and 
Volovik ^ where they considered a moving texture or domain wall in 3 He. 
The sound speed within the wall is different from that outside. And the wall 
needs to be moved relatively slowly with respect to the container for there to 
be Hawking radiation. One complication there is that the order parameter 
varies within the wall in such a way that there is also radiation due to other 
effects. 

The basic experimental set-up seems very simple (see Fig.|5J). There is a 
tube with the fluid in one phase on one side and a second phase on the other 
side. We assume that the sound speed in phase 2 is less than the sound speed 
in phase 1: c\ > ci. If the phase boundary propagates from phase 2 into 
phase 1 at a speed v such that c\ > v > C2, sound from within the phase 2 
region cannot enter the phase 1 region. Alternately, in the rest frame of the 
boundary, the fluid is flowing with velocity v toward the phase 2 side, and the 
Unruh set-up is exactly duplicated. The phase boundary is the location of 
the sonic horizon. Since the fluid is not moving with respect to the container, 
there won't be any instabilities due to surface effects. Although, the phase 
boundary does move with respect to the container and one needs to worry 
if this will lead to instabilities. 
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Fig. 2. A supercooled or superheated system with two phases can be a black hole 
analog if the phase boundary propagates with subsonic speed with respect to Phase 
1 but supersonic speed with respect to Phase 2 (ci > v > c 2 ). In the rest frame of 
the boundary, the fluid is moving to the right at a speed v that is subsonic to the 
left of the boundary and supersonic to its right. 



The Hawking temperature for the propagating phase boundary can be 
calculated in the rest frame of the boundary. Then the metric is of the 
Painleve-Gullstrand-Lemaitre form (Eq. (j2J)) and one can directly use the 
calculations in Refs. 17151 in which Hawking radiation is viewed as the tunnel- 
ing of ingoing particles from inside the black hole to outgoing particles on 
the outside. The advantage of this approach is that it is largely system inde- 
pendent. It does require, however, that the particles be massless both inside 
and outside the black hole. Let us now outline this calculation following 
Ref.|HJ 

In both the exterior and interior regions of the sonic hole (Phases 1 and 
2 in Fig. [2J), the dispersion relation is assumed to be that for a massless 
fermionic particle. In the fluid rest frame: 

E(p) = ±c\p , exterior (6) 

E(p) = ±C2P , interior (7) 
For convenience we will write these equations as: 

E(p) = ±cp (8) 

where c(x) is c\ in the exterior region and C2 in the interior region. The 
transition in c(x) is assumed to be sharp compared to any other length scale 
of interest but c(x) is still assumed to be continuous. 

As for any fermion, the dispersion relation has states with both negative 
and positive energy. The negative energy levels are all occupied and form 
the Dirac sea. 

If we work in the rest frame of the phase boundary, the fluid moves to 
the right with speed v in Fig. |2 Now the velocity of a quasiparticle, defined 
as V p S, is shifted by the fluid velocity. So the dispersion relation becomes: 



E(p) = (±c + v)p 



(9) 
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Fig. 3. The dispersion relation outside (left diagram) and inside (right diagram) 
the dumbhole. Due to the supersonic flow inside the dumbhole the occupied levels 
in the Dirac sea (dashed lines) emerge and get positive energy. Now particles in the 
emerged states can tunnel to an unoccupied state in the exterior region 13. 



The usual situation is when v < C2 < c\. Then the filled levels (Dirac 
sea) of both branches of the dispersion relation have negative energy. What 
happens if v > C2? Then one branch of the Dirac sea in the interior region 
has positive energy as illustrated in Fig- El This is what makes it possible for 
particles of energy E in the interior vacuum to emerge out of the dumbhole. 
As shown in Fig. |HJ it is the particles that have velocity directed to the right 
(further into the interior of the dumbhole), that can tunnel out onto the 
branch for the left-moving particles in the exterior region and then escape 
out to the asymptotic region. 

The energy is related to the momentum of the particle via E = (±c+v)p. 
We are interested in the particle that is outgoing in the exterior region, 
therefore we choose the — sign. 

The Hawking temperature is calculated by finding the tunneling rate, 
which can be obtained from the imaginary part of the action: 

/+oo r+oo 
p(x)dx = Im / —dx (10) 
-oo J-oo V - C(X) 

The function v — c{x) has a zero at the location of the phase boundary 
(x = 0). So we need to go over to the complex plane, deform the contour 
of integration around the pole at the origin, and then evaluate the integral 
using the residue theorem. This gives: 



(11) 
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where |c'| is the magnitude of the derivative of the sound speed at the location 
of the horizon. The rate at which particles can tunnel to the exterior region 
is given by exp(2 lmS) = exp(E /Th) where Th = \c'\/2ir is the Hawking 
temperature. 

Then the formula for the Hawking temperature is the same as in Eq. Q). 
The Hawking temperature is estimated to be: 



where, £ is the thickness of the phase boundary. 

Note that v enters the estimate since it determines the location of the 
horizon at which gradient of c is evaluated. However, we will assume that 
T s h is (roughly) independent of v as long as it satisfies c\ > v > ci . 

The hardest problem, of course, is to find an experimental realization of 
a dumbhole. Since the Hawking temperature is very low, the system must 
either be a fluid or a solid. Only a few fluids are known at this temperature. 
Solids may be used too if there is a suitable melting transition. Another 
possibility is to use Bose-Einstein condensates that have been discussed in 
other related contexts in Refs. I9|10pill 

One particularly intriguing system is superfluid 3 He since the AB phase 
boundary has been studied quite extensively both experimentally and the- 
oretically. Indeed the AB interface has been made to oscillate by applying 
suitable magnetic fields, and the Lancaster group has also measured the ra- 
diation from the oscillating phase boundary^. The experimental setup is 
essentially that shown in Fig. [21 except that the container is vertical. The 
fluid is kept at 150 fiK and bar pressure. At such low temperatures there 
are few excitations present and the system is essentially in its vacuum state. 
A non-uniform magnetic field is applied along the vertical such that it is 
stronger in the lower part of the container than in the upper part. This 
causes the lower part to be in the A-phase and the upper part to be in the 
B-phase. The AB phase boundary is at the location where the magnetic 
field attains a critical value, Bab- Next, a small oscillating magnetic field is 
applied along the vertical so that the total magnetic field is: 



The location of the AB phase boundary in the absence of dissipation can be 
be found by setting B(t, z) = Bab- Due to the oscillating component of the 
magnetic field, the equilibrium location of the interface is time-dependent: 
zq = asin(o;t). The amplitude of oscillation is given by: a = Bac/^Bq. 
(The gradient of Bq is roughly constant in the region where the phase bound- 
ary oscillates.) The frequency of oscillation of the phase boundary can be 




(12) 




(13) 
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Fig. 4. The data showing the dissipation from the moving AB interface as a func- 
tion of the driving frequency as given in Ref. 1 1 121 See text for further explanation. 



changed by changing the frequency, v = uj/2tt, of the applied oscillating 
field. The quasiparticle radiation is detected by a vibrating wire resonator 
placed at the upper end of the container. The data for the radiated power 
versus oscillation frequency is shown in Fig. |1J 

The behavior of the experimental curves in Fig. i3] are explained phe- 
nomenologically as follows. Assuming a linear restoring force and dissipative 
motion, the equation of motion for z(t), the position of the interface, is: 

-k(z-z (t))- 7 ^ = (14) 

where zo(t) = asin(wt) (a = Bac/^JB) is the equilibrium position of the 
interface. Then, 

dz 

v ( t ) = — = v cos(ujt-(f>)+O(e- Kt ) (15) 

where 

k aKu u 

« = -, VQ = — r ===, tan^ = - (16) 

Therefore v o oc u for uj « k and vq is independent of u for uj » k. 
The average power dissipated over an oscillation time period is 

P = l{v 2 ) (17) 



and this grows quadratically with uj for u < k and is constant for u > k. 
This is consistent with the log-log plot of the data for v < 1 Hz in Fig. ^ for 
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the larger value of Bac, an d over the whole range of v for the smaller value 
of Bag provided we treat 7 as a fitting parameter that depends on Bac- 
There are three unexplained features of the data as discussed in Ref. IT2l 

1. Theoretical estimates of 7 are small by several orders of magnitude^. 

2. For the smaller value of Bac, the required value of 7 is about 5 times 
smaller than what is needed to fit the data for larger Bac- If 7 is 
calculated from the Bac = 0.643 mT data, and used to predict the 
Bac = 0.214 mT curve, one gets the lower solid curve shown in Fig. 0] 
and this does not even come close to matching the data. 

3. The data clearly shows an unexpected increase in dissipation at fre- 
quencies v > 1 Hz when Bac = 0.643 mT. 

Can Hawking radiation play any role in the increased dissipation seen 
for v > 1 Hz when Bac = 0.643 mT? The data does seem to indicate a new 
source of dissipation at high frequencies. The velocity of the interface at 
v ~ 1 Hz is about 1 cm/s. This also happens to be close to the quasiparticle 
velocity in the A phase in the direction orthogonal to the interface (ca = 3 
cm/s for certain excitations as explained below). It is also less than the 
quasiparticle velocity in the B phase (cb = 55 m/s). Furthermore, the 
power emitted in Hawking radiation using the temperature in Eq. (|12j) is of 
order 1 pW and is in the range seen. Hawking radiation is only expected 
during the part of the oscillation when the velocity of the interface lies in 
the suitable range: cb > v > ca- As this duration increases, we expect the 
amount of Hawking radiation to grow proportionally. 

At v ~ 1 Hz in the experiment, the simple model for the motion of the 
interface must break down. This could happen in two ways. The restoring 
force might get non-linear corrections, or the dissipation parameter 7 could 
get some velocity dependent corrections. Purely for illustrative purposes, let 
us assume that the form of the velocity (Eq. (|15|)) continues to hold even at 
large to but with a different amplitude: 

v = au) cos(a;t — <fi) (18) 

The amplitude a is assumed to be an u independent parameter in the fre- 
quency range of interest. The exact form of v is not crucial for us. For 
example, the power of u> in the amplitude could be different from 1. 

Let us now assume that Hawking radiation at temperature T s h (Eq. I)12JI ) 
is emitted when v = z > ca- Then, for j "light" species of radiation, this 
gives: 

Ps-i^A^y^Al^^) (19) 
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where St/r is the fraction of time for which v > ca- Note that the formula 
holds only for radiation for which the mass is less than the Hawking tem- 
perature. Since T s h ~ 3mK > Ab ~ 1.7mK this is self-consistent. (If the 
Hawking temperature is less than the mass of the particles, the radiation is 
exponentially suppressed.) Let us write 



D -i 
fjj = — cos 

7T 



(20) 



where 



and d is the cell diameter (4.3 mm in the experiment). The value of Pq is in 
the observed ranee if £ ~ 300 A. 

Can we estimate the critical frequency z/*? In our illustrative model 



ca 
2na 



(22) 



If a = B AC /VB \ 

interface then the values of are shown in the following 

Table. 
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B A c mT 


Hz 


2.00 


0.643 


14.8 


1.00 


0.643 


7.4 


0.53 


0.643 


3.9 


1.00 


0.214 


22.2 



Clearly the values of z/*, the critical value of the frequency at which Hawk- 
ing radiation starts, do not agree with data. The experiment shows anoma- 
lous radiation starting at some critical frequency, whereas Hawking radiation 
would start at some critical velocity. On the other hand, there is no reason 
to adopt a = Bac/^Bq except that this holds at low frequencies. The 
Hawking radiation explanation can only work if a is independent of VBq at 
high frequencies. However, as we will now see, the interpretation in terms 
of Hawking radiation has additional theoretical issues. 

The main difficulty with the interpretation is that, in our earlier esti- 
mate, we have assumed that the sound quanta are massless on either side of 
the interface, whereas this is not true in the AB system. The quasiparticles 
on the A-phase side of the interface indeed have a linear dispersion relation 
very close to the nodal point. This can be seen from the full dispersion 
relation in the rest frame of the fluid ^ 

1/2 



E(p) 



4( 



ip 



PfY + 



A 2 
Pf 



(23) 
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where vf ~ 55 m/s is the Fermi velocity and pf = m*VF is the Fermi mo- 
mentum with m* ~ 3m# e _3, and A^ ~ 2.02fcsT c where T c is the transition 
temperature from normal to superfluid phased At zero pressure and strong 
magnetic field T c ~ 1 mK. It is easy to check that E(p) vanishes if p = ±pfI- 
Hence the dispersion relation has a node. Now consider excitations near the 
node, p = ±pfI + Sp, with 



1M« Aa 

PF 



i(r 



PFVF 



(24) 



and with momenta perpendicular to I: 5p l 
dispersion relation is linear, 

E = ±c A \5p\ 



. For these excitations the 



(25) 



with ca = Aa/pf ~ 3 cm/s. Therefore we can expect a dumbhole to form 
whenever the interface velocity exceeds 3 cm/s. 

On the B-phase side, the dispersion relation in the rest frame of the 
fluid is: 

E(p) = ±\{e(p) - + A%4] 
L Pf. 



(26) 



where e(p) 
T, 



= ±[{*(P) 

e(pF), and the B-phase gap ~ 1.76/cbTc with 



p 2 /2m*, fi 

, — 0.93 mK at zero pressure and magnetic field. Since E only depends 
2 and it does not vanish at p 2 = pp, there is no node in the dispersion 



on p 

relation. In fact, for p 



\p\ ~ pf we have: 



E(p) 



V F {P-PF 



) 2 + A 2 



B 



1/2 



where we also assume 



As 



P-PF 
PF 



< 1 



(27) 



(28) 



-6 



Note that (A B //x) 2 ~ 10 

Since the B-phase dispersion relation is not of the usual relativistic 
form E = ± \Jp 2 c 2 + m 2 c 4 , it is not possible to think of the quasiparticles as 
propagating on a metric with just the usual second derivative kinetic term. 
However the black hole analogy need not break down! There is still a sonic 
horizon from within which quasiparticles moving at 3 cm/s in the A-phase 

1 Note that the I vectors on the AB interface mostly lie in the plane of the interface. 
Hence excitations that cross the interface necessarily have non-vanishing momenta per- 
pendicular to I. 
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Fig. 5. Vacuum fluctuations in the B phase can produce an // pair, say with / 
having positive energy and / having negative energy. / falls into the dumbholc, 
and if its velocity is very low, cannot escape back out into the B-phase region. / 
escapes and forms Hawking radiation. 

region cannot escape. This is an important ingredient for a dumbhole (see 
Fig. EJ). The present situation is more complicated though, because not all 
quasiparticles in the A-phase move at 3 cm/s. Only the quasiparticles close 
to the nodal point have this velocity. So a negative energy particle that falls 
through the interface and into the A-phase region cannot escape back to the 
B-phase region only if its velocity in the A-phase region is of order 3 cm/s. 

The dispersion relations given above are in the fluid rest frame. In 
the interface rest frame, the fluid is moving to the right with velocity v. 
Quasiparticle velocities in the interface rest frame will be shifted by v with 
respect to the velocities in the fluid rest frame. Since velocity is defined 
by 'VpE, the dispersion relation in the interface rest frame is obtained by 
adding v ■ p to the dispersion relation in the fluid rest frame. As before, 
part of the Dirac sea on the A-phase side of the interface will emerge into 
the positive energy region. Furthermore, the energy of some occupied states 
in the A-phase side is the same as that of unoccupied states in the B-phase 
side. So the basic set-up of the linear dispersion relation calculation given 
above is available. The technical difficulty is that the equation e = E(p, c{x)) 
cannot as simply be inverted to get p(E, c(x)). Hence let us construct a toy 
model for the dispersion relation that captures some of the essential features 
of the full model. 

Consider the one-dimensional dispersion relation: 



where c = vp, A = A# in the B-phase region and c = ca, A = in the A- 
phase. We will take P = pp = Pb in the B-phase and P = in the A-phase, 
but also keep open the possibility that Pb = 0. The functions c(x), A(x) 
and P(x) only vary appreciably within the interface. This dispersion relation 
goes over to Eq. (|25|) in the A-phase (with v = 0) and likewise Eq. I|27|l in 
the B-phase 2 . Note that the dispersion relation is written in the rest frame 

2 Though the linear dispersion in the A-phase is due to the A\ term and not due to the 
kinetic term as in the toy model. 




(29) 
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of the interface in which the fluid is moving to the right with speed v. Also, 
we are looking at the one dimensional problem where the momentum is only 
along the x direction. 

Inverting Eq. we get, 



(v 2 — c 2 



P~P= Y ,J ^ " (30) 



where e = E — vP. The discriminant should be positive in the asymptotic 
regions for the particle to be physical (and not "under the barrier"). This 
imposes a constraint on the velocity: 

^ ^ (31) 

In the A-phase region, A = 0, and the condition is trivially satisfied as seen 
from Eq. pT|). In the B-ph ase, the condition reads: 




2 . 2 
v > v F 



1 _ {E-vPb? 



A 2 
B 



(32) 



Hawking radiation can only occur if: 



E>^ B ^l--j T + vP B (33) 

or else if: 

E < -A B] jl-^- + vP B (34) 

We assume that the condition in Eq. (|31|) is satisfied everywhere, including 
within the interface. If this is not the case, the discussion will be more 
involved. 

In the experiment, v ~ ca ~ 10 _3 vf and so the square root factor is 
essentially 1. Therefore either E > Ab + vPb or else E < —Ab + vPb- Note 
that we are interested in v > Aa/pf = ca and since ~ Ab, we also have 
vpp > Ab- So the thresholds for E are both positive if Pb = Pf (but not 
if, for example, Pb = 0). The first possibility is one that is expected since 
the minimum energy of quasiparticles in the B-phase is Ab in the fluid rest 
frame. The second possibility corresponds to the emergence of the Dirac sea 
in the B-phase when viewed in the rest frame of the interface. Let us discuss 
both possibilities in some more detail. 

At the large values of the energy required by the first possibility (Eq. (|33j)) . 
our toy model starts becoming suspect and a more realistic calculation is 
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called for. However, if we proceed with the toy model under the assump- 
tion that it still gives the correct qualitative behavior, we can calculate the 
Hawking temperature. The calculation of the action follows the linear case 
described earlier. There is a pole at v = c{x) and this gives the only contri- 
bution to the imaginary part of the action. However, the process requires a 
branch change. The particles on the A-phase side are on the branch with the 
— sign in Eq. IJ29JI while, after tunneling, the particles would end up on the 
branch on the B-phase side that has the + sign. Further analysis is required 
to determine if the branch change can take place. 

Now we discuss the second possibility (Eq. (|34j)). In this case, the 
particle from the A-phase may tunnel to a state in the Dirac sea of the B- 
phase. However, the tunneling cannot take place since the corresponding 
state in the B-phase Dirac sea is occupied. Yet the process may still be 
important. Suppose there is a vacuum fluctuation in the B-phase in which 
a particle from the Dirac sea gets excited to the upper branch. Normally 
the particle would fall back into the hole in the Dirac sea within a time 
allowed by the uncertainty principle. However, in the present situation, a 
particle from the A-phase can tunnel into the B-phase and fill up the hole 
before the original particle has had a chance to fall back. Then there is 
no hole left for the particle in the upper branch to return to and it must 
escape as a real particle (see Fig- EJ) - This literally corresponds to the process 
shown in Fig. |S] where vacuum fluctuations in the B-phase create a pair of 
positive and negative energy fermions and the negative energy particle falls 
into the dumbhole. So the second possibility given by Eq. (|3*l)l might indeed 
be relevant. 

The calculation of the Hawking temperature proceeds as before. The 
tunneling rate is calculated by finding the imaginary part of the action ob- 
tained by integrating q(x). Note that the + sign in Eq. (|3Uj) must be chosen 
for there to be a pole in q(x). The tunneling rate must be multiplied by the 
probability of having a hole of energy E on the B-phase side. As long as this 
factor is not an exponential, the Boltzmann factor will be given by the tun- 
neling rate alone. Since the structure of the pole is the same as in the linear 
calculation given above, the Hawking temperature will still be T s h = \c'\/2ir 
in the interface rest frame. Note that since the tunneling occurs between the 
branches with a — sign in Eq. ()29|) . E must lie in the range (for Pg = Pf)'- 

< E <v Pf - A B (35) 

where we are assuming v <C vp. Therefore, in the fluid rest frame, the 
radiated particle can only have energy (E r ) in the interval: 



As < E r < vpf 



(36) 
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Fig. 6. The figure shows the schematics of the dispersion relations in the toy model 
for the A- and B- phases. Dashed curves denote the filled states of the Dirac sea. 
If there is a vacuum fluctuation in the B-phase, a hole can be created in the lower 
branch, which can then be re-filled by a particle tunneling out of the A-phase. The 
particle that had jumped onto the upper branch during the vacuum fluctuation then 
escapes to infinity as Hawking radiation. This process is identical to that shown in 
Fig.E 



This restriction on the range of energy of radiated particles will change if 
branch changing processes can occur. Then one can imagine several other 
processes as well. These should be investigated. 

The above discussion is based on a toy model of the AB system and 
suffers from the danger that it is inaccurate in some essential way. In par- 
ticular, since the last term in Eq. (|2 u fl comes from v ■ p, and the component 
of p along v is taken to be small on the A-phase side, perhaps it would have 
been better to replace the term by v(p — P(x)). Then, with Pg = pp and 
for small p — pp, the B-phase Dirac sea would not emerge into the positive 
energy region and the only possibility for tunneling would be via branch 
changing processes. 

Clearly it is important to do a more thorough calculation and with a 
realistic model of the AB interface. Such an effort would be aided by the early 
calculations on the dissipation from the AB interface the profile of the 
order parameter within the interface and the scattering of quasiparticles 
off the interface To connect with experiment, one would also need a 
handle on the interface motion at high frequencies and on other effects that 
could potentially explain the anomalous radiation. 

Experimentally, from the perspective of studying Hawking radiation, it 
would be desirable to obtain the spectrum of the emitted radiation, and to 
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be able to correlate the emitted power with the dynamics (say, the velocity) 
of the interface. 

If it turns out that Hawking radiation from the AB interface of 3 He is 
highly suppressed (or absent), we will be left with the challenge of finding a 
quantum black hole analog. This is an exciting quest. It is probably also our 
only hope for experimentally testing current ideas on quantum black holes 
and other cosmological problems. 
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